A new solution to the solar neutrino deficit 
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We reexamine the transition magnetic moment solution to the solar neutrino problem. We argue 
that the absence of large time variations in the Super-Kamiokande rate provides strong evidence 
against spin-flavor flip in the solar convective zone. Spin-flavor flip could, however, occur in the 
primordial magnetic field in the radiative zone. We compute the longest lived toroidal mode for this 
field and show that spin-flavor flip in the radiative zone can account for all available solar data. 
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I. INTRODUCTION 

In recent years, the progress in the field of solar neu- 
trino physics has been remarkable. This progress culmi- 
nated when the SNO collaboration conclusively demon- 
strated the presence of a non-electron neutrino compo- 
nent in the flux of solar neutrinos [1,2]. 

The SNO experiment measured the flux of the 8 B solar 
neutrinos using three different reactions: (i) the charged 
current (CC) process v e + d — ► p + p + e~ , (ii) the elas- 
tic scattering (ES) process v x {v x ) + er — > v x {v x ) + e~ , 
and (iii) the neutral current (NC) process v x (v x ) + d — > 
p + n-\-v x (p x ). A high statistics ES measurement, in the 
same energy range, had also been previously carried out 
by the Super-Kamiokande (SK) experiment [3]. While 
the CC process selects only electron neutrinos, the other 
two processes are sensitive to neutrinos as well as an- 
tineutrinos of all active flavors. 

The results of the measurements are as follows [2,4]:* 



i?cC Q _ i.7 6 ± o.ll x 10 6 cm _2 s _1 , 
Rsno = 2 - 39 ± °- 27 x 10 6 cm- V 1 , 
= 2.35 ± 0.08 x 10 6 cm _2 s -1 , 



tfNCi = 5 09 ± 0,63 x 10 6 cm _2 s -1 , 
i?NC2 = 6.42 ± 1.67 x 10 6 cm _2 s -1 . 



(1) 
(2) 
(3) 
(4) 
(5) 



They are impossible to reconcile assuming the incident 
flux consists entirely of electron neutrinos [1,5,2], and 
hence some other particles must be contributing to the 
ES and NC event rates. Chronologically, the first detec- 
tor to observe events caused by these other particles was 
SK, but it was not possible to separate these events from 
the f e -induced events until the SNO CC data became 
available. One can therefore say that the SNO measure- 
ment has effectively reclassified SuperKamiokande as an 
appearance experiment. 

The interpretation of the SNO results, however, is not 
unambiguous. While the "excess" NC and ES events are 



*J? SNO was found assuming undistorted B neutrino energy 
spectrum, while -Rsno was f° un( l when this assumption was 
relaxed. 



most often ascribed to muon (and tau) neutrinos, they 
can be also caused by muon (and tau) tmimeutrinos. 
The antineutrinos in question must be of the muon or 
tau types, since electron antineutrinos would be easily 
identified through the reaction v e + p — ► n + e + [6]. 
Both v^ r and P M;T scatter on electrons and deuterium 
nuclei through their neutral current (NC) interactions, 
with similar cross sections. At E v ~ 10 McV, 



da(D^e-) 



dT, 



d/j(i 



dT e . 



0.8, (6) 



where the integral over the electron recoil energy T e 
ranges from T min = 5 MeV to T max = E v /{1 + m e /2Ej) 
MeV. Similarly, using the latest available calculations of 
the v x (p x ) + d —>■ p + n + v x {p x ) cross sections [7] and 
integrating over the undistorted 8 B reaction spectrum, 
f{E v ), [8], one finds 



j dE v f(E v )a{D^ r d)/ J dE v f{E v )a{y^ T d) ~ 0.95. 



(7) 



Assuming the v e — > v^ T conversion mode, Eqs. (1-4) 
imply that the electron neutrino survival probability is 
P ee ^ 0.35 ± 0.04 {Pee ^ 0.27 ± 0.06 if Eqs. (1-3,5) are 
used). If instead the v e — > v^,t mode is assumed, the 
same data imply a somewhat higher flux of all active 
neutrinos, namely, <E>tot = 5.25 ± 0.66 cm -2 s _1 if one 
uses Eqs. (1-4) or $ tot = 6.64± 1.75 cm" 2 s _1 if one uses 
Eqs. (1-3,5). The corresponding survival probabilities are 
Pee ^ 0.34 ± 0.05 or P ee ^ 0.27 ± 0.07. Thus, both con- 
version modes are capable of accounting for the measured 
event rates. Moreover, both yield the total flux that is 
consistent with the Standard Solar Model (SSM), which 
predicts 5.05±J;g x 10 6 cm" 2 s^ 1 [9] or 5.93 ±0.89 x 10 6 
cm~ 2 s _1 [10], depending on the input value of the S17 
parameter used. Therefore, in interpreting the experi- 
mental data both possibilities should be considered. 



The 



conversion is predicted in the flavor 



oscillation scenario [11], while v e — > v„ T is predicted in 
the neutrino spin-flavor flip (SFF) scenario [12-16]. The 
flavor oscillation scenario has been extensively studied 
over the years. The goal of this paper is to reexamine 
the SFF scenario and to point out the existence of a new 
solution. 
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In the SFF scenario, the neutrino is postulated to cou- 
ple to the electromagnetic field through a magnetic mo- 
ment, //,, interaction. A neutrino propagating in the solar 
magnetic field undergoes a spin rotation, and, if it is a 
Majorana fermion, becomes an antincutrino, in the pro- 
cess also changing its flavor. 

Majorana fermions are most naturally described using 
the two-component Weyl spinor notation. The coupling 
in question is given by a Lagrangian term 



C 



EM 



(8) 



where (<r* 



P = 



- T af3 v 



aa a a 



dtp p 



)/2, 



= (1,<t) 



(1,— ct). \i and v are Lorentz indices, a and (3 
are Weyl spinor indices, and a and b are the flavor in- 
dices. Since the spinors anticommute, the a = b terms in 
Eq. (8) vanish identically, and hence Majorana neutrinos 
cannot have magnetic moments [17]. They can, however, 
have transition (a b) moments which lead to a simul- 
taneous spin and flavor change. This provides a natural 
scenario for v e — > v^^ transitions without also produc- 
ing v e states. Flavor mixing, if present, must be small to 
avoid a significant secondary v^^ — ► v e conversion. Re- 
garding the naturalness of this assumption we note that 
out of the three mixing angles in the leptonic sector the 
one measured with atmospheric neutrinos is large and the 
one measured by reactor experiments is small, so that a 
priori there is no natural value for the third angle. 

In general, the evolution of the neutrino state is gov- 
erned by a 6 x 6 Hamiltonian matrix. For simplicity, we 
consider a 4 x 4 case. In the basis (v e , v^, P e , 9^), the 
Hamiltonian is [15] 



H 



H v (B x - iB y )M^ 

(B x + iB y )M H D 



(9) 



where B x ^ y are the transverse components of the mag- 
netic field and the 2x2 submatrices are given by 



M = 
H v = 

H B = 









-Acos26» + A e A sin 26 
A sin 26 A cos 20 + A 



-A cos 26 - A e 
A sin 20 



A sin 26 
A cos 26 - A u 



(10) 

(11) 
(12) 



Here A = Am 2 /AE V , A e = V2G F (n e - n„/2) and 
= v / 2Gf(— n n /2). Am 2 is the neutrino mass-squared 
splitting, E v is its energy, and n e and n n are the electron 
and neutron number densities. 

As will be described later, to obtain a significant 
v e — * Vft.r conversion in the Sun requires fi e ^ <~ (10~ 12 — 
10 _10 )^b (fiB = e/2m e is the Bohr magneton). The 
transition moment values in this range are below the di- 
rect laboratory bounds derived from the analyses of the 
v + e — > v + e [18] and v + e — > D + e [19] scattering 



data. Stronger bounds on the transition moment have 
been claimed from the analysis of various astrophysical 
processes, most notably from the study of the red giant 
populations in globular clusters [20]. The idea is that 
a sufficiently large transition moment would provide an 
additional cooling mechanism and change the red giant 
core mass at helium flash beyond what is observation- 
ally allowed. The bound stated in [20] is < 3 x 10~ 12 ^b. 
Recently, there has been some disagreement between the 
stellar models and also a systematic shift to longer dis- 
tances for clusters [21], and it would be very interesting 
to see an updated bound. 

The compilation of various bounds on the size of the 
transition moment and the corresponding references can 
be found in [22]. 



II. MAGNETIC FIELD IN THE CONVECTIVE 
ZONE 

The neutrino spin-flavor flip can occur either in the 
convective zone (CZ) (0.71i?o < r < Rq) or in the ra- 
diative zone (RZ) deeper in the solar interior. The two 
possibilities are physically quite different and require sep- 
arate treatments. First, we consider the CZ case. 

While the magnetic field in the inner part of the CZ 
cannot yet be measured directly, its basic properties can 
be inferred from the field at the solar surface, and it 
is widely believed that large scale field structures exist 
at some depth in the CZ [23]. Recent magnetohydrody- 
namic dynamo models argue that the field generation oc- 
curs in the shear layer near the bottom of the CZ. These 
models predict field values as large as 100 kG [23], which 
significantly exceeds the turbulent equipartition value of 
B - pVe^VV 2 / 3 - 10 4 G t. The helioseismological 
data provide an upper limit of 300 kG on the magnitude 
of this field [24]. 

The field strength necessary for significant v e — > v^. T 
conversion can be estimated as fi e ^BAl <~ 1, where Al is 
the thickness of the magnetic field. One finds 



(^/hb){B/1 kG)(A//0.1i? Q ) - 5 x 10 



-10 



(13) 



so that for [i ell <~ 10 _11 ^s and Al ~ 0.1i? Q the value 
B <~ 50 kG is required. Although the surface field, which 
reaches several kG in sunspots, is too weak to affect the 
evolution of the solar neutrinos, the field near the bot- 
tom of the CZ might have the right strength (provided 
fj, elJl ~ 10~ n /is). This observation, and the natural ex- 
pectation that the field in the interior of the CZ changes 



^This estimate is obtained by estimating the convective en- 
ergy flux as Lq ~ 4irr 2 vpv 2 /2, where v is the turbulent ve- 
locity, and assuming equipartition pv 2 /2 ~ B 2 /Stt. 
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with the solar cycle, were the reasons the convective zone 
SFF were originally suggested as an explanation for an 
apparent anticorrelation between the Homestake event 
rate and the number of sunspots [13]. 

Over time, however, a paradigm shift has occurred. 
The idea of using this mechanism to explain variations 
of the neutrino rates was abandoned, especially in view 
of the SK data, which showed no seasonal or yearly vari- 
ation, except for a hint of the expected 1/r 2 flux mod- 
ulation. Instead, recent analyses [25] assume that the 
magnetic field in the CZ is somehow time independent, 
an assumption that we find untenable. The presence of 
a strong constant magnetic field inside the CZ would, on 
the one hand, almost certainly be revealed at the surface 
by convective mixing and, on the other hand, cause an 
asymmetry in the solar cycle. The strong toroidal field 
that is believed to exist in the overshoot layer, where con- 
vective mixing is small and the shear is large, should also 
be variable because it is generated from a variable radial 
field. It is therefore implausible that the SFF process in 
the CZ could cause a ~ 70% depletion of the 8 B neutrino 
flux, while producing no observable variations at SK dur- 
ing more than four years. This time span covers about 
half of the solar cycle, with an average sunspot number 
- 10 in 1996, and - 120 in 2000 [26]. 

In addition to the year-to-year variations, one also ex- 
pects characteristic semiannual variations [13], for the 
following reason. The field on the surface of the Sun 
changes direction across the solar equator and vanishes 
at the equator (where indeed no sunspots are observed). 
The transition region, where the field is small, has the 
size ~ 7 x 10 9 cm [13,27]. The orbit of the Earth 
makes a 7° angle with the plane of the solar equator 
and, since 95% of the 8 B neutrinos are produced in 
the region r < 0.09i? Q [9], a significant fraction (88%) 
of them passes through the transition region when the 
Earth crosses the plane of the solar equator. On the 
other hand, in March and September, when the Earth 
is farthest from the solar equatorial plane, this fraction 
is as low as 1%. Thus, large semiannual variations are 
expected at Super-Kamiokande and are not observed [3] . 

We therefore regard the absence of variations in the 
SK data as strong evidence against SFF in the CZ. 

III. RELIC FIELD IN THE RZ 

We now turn to the possibility that the SFF process 
could occur in the radiative zone of the Sun. Unlike the 
convective zone, the radiative zone is not continuously 
mixed and rotates as a solid body. Therefore, it could in 
principle support a "frozen" magnetic field configuration. 

The properties of magnetic fields in the solar interior 
were investigated by Cowling [28] , and by Bahcall and Ul- 
rich [29]. These authors consider poloidal magnetic fields. 
We have redone the calculation [30] using a modern solar 




FIG. 1. The radial dependence of the field strength for the 
lowest toroidal mode in the RZ of the BP2000 SSM. 



model (BP2000), and, more importantly, assuming that 
the field is toroidal. A strong poloidal field would pene- 
trate the convective zone and would either dissipate over 
the solar lifetime, or be observable. We also assumed 
axisymmetry, which is physically plausible. Large asym- 
metry would lead to observable 27-day variations of the 
neutrino rate due to solar rotation. 

For an axisymmetric toroidal field, the magnetic dif- 
fusion equation (describing Ohmic decay) was solved 
numerically. The eigenmodes are of the form B = 
e -t/n,„^ n ( r )pl( cos e) ) w here / = 1,2,.., n = 1,2,.., 
P/^cosO) is the associated Legendre polynomial and 
F (r) is determined numerically. The I = 1, n = 1 mode 
has the longest lifetime, t\ i — 24 Gyr; its radial pro- 
file is shown in Fig. 1. For higher modes, T2,i = 13 Gyr, 
Ti,2 = 10 Gyr, ... For comparison, the solar age is thought 
to be 4.6 Gyr [9]. Therefore, a toroidal field of complex 
spatial structure can exist in the RZ of the Sun. 

Observations provide upper bounds on the allowed 
strength of this field [30]. Neutrino fluxes and helio- 
seismological measurements of the sound speed both 
give B < 10 2 MG. Possible residual (non-rotational) 
oblateness and helioseismological measurements of the 
frequency splittings provide stronger bounds, B < 7 MG 
[30]. 

We next consider the neutrino evolution in the pres- 
ence of this relic field and show that with appropriate 
parameter choice there exists a fit to all available solar 
neutrino data. 

For simplicity, we assume no flavor mixing and that 
only the longest living eigenmode is present. This is suf- 
ficient to demonstrate the existence of the solution. The 
neutrino evolution equations then reduce to a pair of cou- 
pled equations, 

idiip e = (A e (r) - A)ip e + n el ^Bx(l)ipp,, 

idii/>n = n e ^B ± (l)^ e + (A - A^r))^. (14) 

Here I is the distance along the neutrino trajectory and 
B± is the transverse component of the magnetic field. 
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2. Electron neutrino survival probability P ee as a 



function of B max and Am /E v , for fj, efi = 10 11 /is- 



There are two free parameters, Am 2 and the product of 
the neutrino magnetic moment fj, efJ , and the normaliza- 
tion of the magnetic field B max . For dcfinitivencss, we 
set /z e/J — 10~ n [Ib- Later we will address the possibility 
of having different values of ix ell . 

Fig. 2 shows the survival probability P ee for electron 
neutrinos produced in the center of the Sun, as a function 
of -B max and Am 2 / E u . The probability was obtained by 
solving the evolution equations (14) numerically on a grid 
of points. 

The behavior of P ee can be given a simple analytical 
interpretation. The situation is completely analogous to 
the classical case of the small angle MSW effect [15,14]. 
The SFF process takes place in a thin layer around r res 
given by the condition A e (r ICS ) + A^(r vcs ) — 2 A. For 
AmtjEv > 10~ 5 eV 2 /MeV, the resonance never oc- 
curs and therefore P ee — > 1. For Am 2 / E v < 4 x 1CT 8 
eV 2 /MeV, the resonance occurs close to the edge of the 
RZ where the field is small and hence once again P ee — > 1 
(or outside in the CZ, in which case P ee varies with time). 
In the intermediate region, the conversion efficiency de- 
pends on the value of fi epi B±^ at ? 



; r rcs : for fi ell B± 
> 0.7 x 10~ n MGx/ib the conversion is adiabatic and 

Pee 0. 

Notice that for this description it is crucial that every- 
where except in the thin resonance layer the off-diagonal 
term fi e p,B± is much smaller than the diagonal splitting 
(A e (r)+A fl (r))/2 — A. This allows one to ignore the fact 
that the off-diagonal coupling varies along the neutrino 
trajectory and use a small angle MSW formula for the 
level crossing probability 



Pr 



exp 



2(^B ± (l)) 2 



\d(A e + A^/dl\ 



(15) 



The complications related to the large angle MSW ef- 
fect [31] are likewise avoided. The neutrino survival 
probability is given by P ee = (1 + (1 — 2P C ) cos20 )/2, 
where cos26> = (2A - A e - A^)[(2A - A e - A^) 2 + 
{2[i elL B i) 2 ] -1 / 2 is computed at the production point. 
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FIG. 3. Regions in the (Am 2 , _B max ) parameter space 
where the rates of various solar neutrino experiments are re- 
produced. 



We performed a fit to all available solar neutrino data. 
In our computations, we used the fluxes of the pp, 7 Be, 
8 B, 13 N, 15 0, and pep neutrinos, as predicted by the 
BP2000 solar model [9]. The data to be fitted include 
the average rates of the GALLEX/GNO [32] and SAGE 
[33] experiments, 

Rqa = 74.7 ± 5.1 SNU, (16) 

the rate of the Homestake experiment [34] , 

Rchl = 2.56 ± 0.23 SNU, (17) 

the CC and NC rates of the SNO experiment as given in 
Eqs. (1,4), and the rate and the energy spectrum of the 
SK experiment [3]. 

For a given production point in the solar core, the 
neutrino survival probability is determined according to 
Eqs. (14). To obtain the total flux of the neutrinos of 
a given energy, one needs to integrate over the neutrino 
production regions. The problem is technically quite non- 
trivial, because of the 3-dimensional structure of the field 
and because only the field component transverse to the 
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2 4 6 8 10 12 14 
E e (MeV) 

FIG. 4. Survival probabilities for an electron neutrino born 
in the center of the Sun for three values of B max : 0.45 MG, 
0.5 MG, and 0.55 MG (top to bottom). Am 2 = 2.5 x 10" 6 
eV 2 was chosen. 



neutrino trajectory enters the evolution equations. * The 
correct integration is absolutely essential for computing 
the converted fraction of the low energy pp neutrinos, 
which have a broad production region and for which the 
resonance occurs close to, and even inside, that region. 

Fig. 3 shows regions of the parameter space where 
the event rates predicted for various experiments agree 
with the corresponding measured rates. The dashed 
lines delineate the bands where agreement is reached 
with a particular experiment. The intersection region, 
where a good fit to all rates is obtained, is character- 
ized by 0.45 MG< B max < 0.55 MG and 1.5 x 10~ 6 
cV 2 < Am 2 < 3 x 10~ 6 cV 2 . It represents our new solu- 
tion to the solar neutrino deficit. 

As an illustration, we list below our predictions for the 
point Am 2 = 2.5 x 10~ 6 eV 2 and B max = 0.5 MG. We 
find that for this point the 8 B flux is suppressed by a 
factor of 0.34, the pp flux - by a factor of 0.86, and the 
7 Be flux - by a factor of 0.36. The 7 Be suppression factor 
implies the effective Borexino flux suppression of 0.44, 
which is below 0.65^$ 12 expected for the LMA solution 
[10]. The predicted rates are 80.6 ± 2.6 SNU for the Ga 
experiments and 2.69 ± 0.31 SNU for the CI experiment. 
The corresponding \ 2 for the rates fit is 3.4 for 4 d.o.f. 

In Fig. 4 we show the survival probability for a neu- 
trino produced in the center of the Sun as a function of 
the neutrino energy, taking Am 2 = 2.5 x 10 -6 eV 2 . The 
three curves correspond to three different normalizations 
of the magnetic field eigenmode, B max = 0.45,0.5,0.55 
MG. We caution the reader that the flux of the pp neu- 



* Fortunately, in the case of the lowest field eigenmode a 
considerable simplification occurs because of the sin O depen- 
dence of the field on the polar angle. Still, a 2-dimensional 
numerical integration is required. 
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FIG. 5. Predicted SuperKamiokande recoil electron energy 
spectrum for B max = 0.50 MG and Am 2 = 2.5 x 10" 6 eV 2 . 



trinos cannot be easily read off from this plot, since, as 
already mentioned, to compute it one has to appropri- 
ately integrate over the production region. 

The survival probability is not strictly constant for 
high energy 8 B neutrinos, however, its variation with en- 
ergy is rather small. The recoil electron energy spectrum 
expected at SK is shown in Fig. 5. The value of the \ 2 
is 24 for 19 d. o. f. The spectrum becomes flatter as one 
increases Am 2 . 

It is worth stressing that the value of the field that gives 
a good fit trivially scales with the value of the transition 
moment according to B' ma-X = -Bmax X 10 -1 Vs//i eAI . In 
particular, if one chooses fi elJ , = 3 x 10~ 12 ^b in order 
to satisfy the astrophysical bound given in [20], one can 
choose B' max ~ 1.5 — 2 MG, which is still perfectly ac- 
ceptable, as shown in [30] . No property of the field, such 
as its lifetime, etc., is affected by such a rescaling. 

Our scenario may be convincingly tested (and perhaps 
excluded) by the KamLAND experiment, should it es- 
tablish flavor oscillations with the LMA parameters, or 
by the Borexino experiment, should it observe large day- 
night variations, as predicted for the LOW solution (SFF 
predicts no variation). 

There are three potential sources of time variations in 
our scenario. First, the effect of the time varying CZ field 
can be shown to be negligibly small, assuming Bcz ^ 100 
kG and Buz and Am 2 in the range of the fit. Second, 
if the field is not axially symmetric, one expects 27-day 
variations. Third, if, in addition to I — 1 modes, the 
I = 2 modes (P 2 1 (cosO) oc sin 20) are also present, the 
flux may exhibit annual variations, for the same reason 
that semiannual variations are expected for the CZ SFF. 
These variations have extrema in March and September, 
and hence are distinguishable from the variations due to 
the eccentricity of the Earth's orbit. Calculations show 
that for the points in the allowed region a relative change 
in the SK rate is 5R/R ~ 2SB/B. Hence, even if the I = 2 
and I = 1 modes have comparable amplitudes, the rate 
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variation is expected to be of the order of a few percent, in 
contrast to SFF in the CZ, in which case it is expected to 
be much larger. The current data provide a weak upper 
bound on the strength of the I = 2 modes. Both the 27- 
day and the annual variations described above constitute 
smoking gun signatures of the RZ SFF mechanism. 

In summary, our scenario provides a fit to all avail- 
able solar neutrino data, is compatible with all existing 
constraints, and may be tested in the very near future. 
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